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Within hydrodynamics we study the effects of the initial 
spatial anisotropy in non-central heavy-ion collisions on the 
momentum distributions of the emitted hadrons. We show 
that the elliptic flow measured at midrapidity in 158 A GeV/c 
Pb-|-Pb collisions can be quantitatively reproduced by hydro- 
dynamic expansion, indicating early thermalization in the col- 
lision. We predict the excitation functions of the 2"^* and 4"^ 
harmonic flow coefficients from AGS to LHC energies and 
discuss their sensitivity to the quark-hadron phase transition. 

The recent observation of transverse collective flow 
phenomena in non-central heavy-ion collisions at ultrarel- 
ativistic beam energies [^]-Q has led to renewed intense 
theoretical interest in this topic (see Q for a review). 
Collective flow is the consequence of pressure in the sys- 
tem and thereby provides access to the equation of state 
of the hot and dense matter ( "fireball" ) formed in the re- 
action zone. This access is indirect since the flow in the 
final state represents a time integral over the pressure 
history of the fireball. Sorge |q] has argued that different 
types of transverse flow (radial, directed, elliptic, see M) 
show different sensitivities to the early and late stages 
of the collision such that a combination of flow observ- 
ables may allow for a more differential investigation of 
the equation of state. In particular, he pointed out that 
the elliptic flow (which develops in non-central collisions 
predominantly at midrapidity and manifests itself as an 
elliptic deformation of the hadronic momentum distribu- 
tions around the beam axis pUq]) is a signature for the 
early stage of the collision: its driving force is the spatial 
eccentricity of the dense nuclear overlap region which, 
if thermalized quickly enough, leads to an anisotropy of 
the pressure gradients which cause the expansion. Since 
the developing anisotropic ffow reduces the eccentricity 
of the fireball, it acts against its own cause and thus shuts 
itself off after some time. Radial flow, on the other hand, 
responds to the absolute magnitude of the pressure gra- 
dients and not only to their anisotropy; it therefore exists 
also in central collisions, and in non-central collisions it 
continues to grow even after the initial elliptic spatial 
deformation of the fireball has disappeared. 

A phase transition from a hadron gas to a color- 
deconfined quark-gluon plasma causes a softening of the 
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equation of state: as the temperature crosses the crit- 
ical value for the phase transition, the energy and en- 
tropy densities increase rapidly while the pressure rises 
slowly. The resulting small ratio oi p/e at the upper end 
of the transition region ("the softest point" Q) weak- 
ens the build-up of fiow as the system passes through 
it. Shuryak [Q and van Hove ||ll| therefore suggested 
that a plot of the mean transverse momentum against 
the central multiplicity density should show a plateau. 
Later hydrodynamic calculations did not confirm the ex- 
istence of a plateau, showing only a slight fiattening of an 
otherwise strictly monotonic curve jlj]. While the accel- 
eration of the matter is weak in the transition region, the 
system also takes a long time to cross it, thereby allowing 
for the fiow to build up over a longer time. This consid- 
erably reduces the sensitivity of the final radial flow to 
the existence of a soft region in the equation of state. 

Recently van Hove's idea was revived in connection 
with elliptic fiow [^ 14|. Sorge [|lj] used a modified 
version of RQMD which allows to simulate an equation 
of state with a "softest point" and found that the re- 
sponse of the elliptic fiow to the initial fireball eccen- 
tricity was weakened for initial conditions in the phase 
transition region. Using a hydrodynamic model, Teaney 
and Shuryak fla] argued that the existence of the phase 
transition should, at higher energies, also lead to other 
dramatic effects in the transverse expansion pattern of 
non-central collisions, in particular to the formation of 
two well-separated shells moving into the reaction plane. 
In the present Letter we follow up on these ideas, trying 
to understand in more detail the transverse dynamics in 
non-central collisions and what experimental data can 
tell us about it. We use a similar hydrodynamic ap- 
proach as in ll^,P, adjust its free parameters to data 
from central Pb-|-Pb collisions at the SPS, demonstrate 
that it correctly reproduces the measured elliptic flow of 
pious and protons at midrapidity |§,0, and then use it 
to make predictions at other beam energies. In particu- 
lar we discuss the sensitivity of the excitation functions 
of V2 and V4, the 2"^^ and 4**^ harmonic flow coefficients, 
to the existence of a deconfining phase transition. 

In the hydrodynamic model one assumes that shortly 
after the impact the produced strongly interacting mat- 
ter reaches a state of local thermal equilibrium and sub- 
sequently expands adiabatically. In the conservation laws 
for energy-momentum and baryon number 



9^T'^''(x) = 0, df,fix)^0 
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one can then use the ideal fluid decompositions T^^'^ = 
{e+p)u^u'^ — g^^p^ j^ = nu'^ in terms of the energy den- 
sity e, the pressure p, the (net) baryon number density 
n, and the fluid four- velocity w'^. One thus obtains the 
equations of ideal (non-dissipative) relativistic hydrody- 
namics. An equation of state p(e, n) is needed to close the 
set of equations; its direct connection with the developing 
flow pattern makes hydrodynamics the most appropriate 
framework for an investigation of the equation of state. 

We are here mostly interested in the transverse expan- 
sion dynamics in non-central {b ^ 0) heavy-ion collisions. 
The lack of azimuthal symmetry leads to a non-trivial 
3-|-l dimensional problem, requiring considerable numer- 
ical resources. As noted in g], the complexity of the task 
is significantly reduced if one focusses on the transverse 
plane at midrapidity and assumes that the longitudi- 
nal expansion can be described analytically by Bjorken's 
scaling solution |l^] Vz — z/t. The latter is known to 
correctly reproduce the longitudinal expansion dynamics 
at asymptotically high beam energies, and it works phe- 
nomenologically very well even at SPS and AGS energies 
p9[. This assumption reduces the numerical problem to 
2-1-1 dimensions. While it should be harmless at midra- 
pidity, it forbids to make reliable predictions at forward 
and backward rapidities. Hence we cannot describe the 
rapidity dependence of the transverse flow pattern. 

We investigated three different equations of state: (i) 
an ideal gas of massless particles, P = § (EOS I); (ii) 
a hadron resonance gas including all known resonances 
pO| with masses below 2 GeV and a repulsive mean 
field potential V(n) = ^Kn"^, with K = 0.45 GeVfm^ 
pl| (EOS H; for small n this equation of state can be 
well characterized by the simple relation p — 0.15 e); 
(iii) an equation of state with a first order phase tran- 
sition at Tc{n=0) = 164 MeV, constructed by matching 
EOS H and EOS I using a bag constant B^/* = 230 MeV 
(EOS Q) 0. EOS Q features at n = a latent heat of 
1.15 GeV/fm^: the mixed phase ranges from en = 0.45 
GeV/inx^ to eg = 1.6 GeV/ini^. We show results only 
for the semi-realistic cases EOS H and EOS Q. 

For 6 7^ the initial energy density distribution in the 
transverse plane has an almond shape, characterized by 
an eccentricity tx = (('^Z^^'^i > 0. {x denotes the trans- 
verse direction parallel to the impact parameter b, and 
((. . .)) indicates an energy density weighted spatial aver- 
age.) This results in larger pressure gradients and thus 
in larger flow velocities in x than in y direction. Hence 
the final px-distribution is anisotropic. Its azimuthal an- 
gular dependence can be characterized by (even) Fourier 
coefficients |1^] W2, f4, . . . (at midrapidity the odd ones, in 
particular the "directed flow" vi, vanish by symmetry): 
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(jj = Artanh(pz/i?) is the longitudinal rapidity of the 
particles, and j/cm below denotes the midrapidity point.) 

For each impact parameter b, we parametrize the initial 
transverse energy density e(r) by a Glauber-inspired for- 
mula which assumes that the deposited energy is propor- 
tional to the number of collisions producing wounded nu- 
cleons (for details see Pp3]). The initial baryon number 
density n{r) is taken proportional to e(r). The propor- 
tionality factors (in particular their -^/s-dependence) can- 
not be calculated but must be adjusted to data. For adia- 
batic hydrodynamic expansion there exists a unique rela- 
tion between the initial entropy density (related to cq, hq 
by the equation of state) and the final multiplicity den- 
sity dN/dy. We therefore present excitation functions as 
functions of the (total) pion multiplicity density at midra- 
pidity, where the "energy calibration" dN^^ / dy\y^y^^^{y^) 
will be provided by experiment. 

The final particle distributions dNi / (dypTdpTdip) are 
calculated using the Cooper-Frye formula |2^ with a 
freeze-out hypersurface of constant temperature. We ad- 
just the model parameters, i.e. the initial central energy 
and baryon densities at b=0, cq and no, the equilibration 
time To, and the decoupling temperature Tjoc, by fitting 
p2[ the measured [M negative hadron and proton spec- 
tra from central 158 A GeV Pb-|-Pb collisions. For EOS Q 
we find Tdcc=120 MeV, eo=9.0 GeV/fm^, no=0.95 fm-^, 
and tq^O.S fm/c. The freeze-out temperature and the 
average radial flow resulting from these initial conditions 
agree well with previous studies p5|-p7[|. In calculating 
the negative hadron spectrum we included |^ decays of 
all resonances up to the mass of the A(1232); resonance 
decays are found to reduce the momentum anisotropics 
W2,4 for pious by 10-15%. 

Having adjusted the model parameters in b—0 colli- 
sions, we can calculate the initial density distributions 
also for 6 7^ collisions, using the same Glauber formula. 
The equilibration time tq and decoupling temperature 
Tdcc are left unchanged. We have tested this procedure 
on pT-spectra from non-central Pb-l-Pb(Au) collisions for 
pions Q and protons [jl|j29| and found very good agree- 
ment between data and hydrodynamic simulations up to 
impact parameters of about 10 fm |g2[. We then proceed 
to compute the 2"^^ and 4'^ harmonic coefficients V2 and 
V4 in Eq. (0) as functions of the number of participat- 
ing nucleons, A^part (or, equivalently, as functions of the 
impact parameter b). The results, for EOS Q and the 
model parameters given above, are shown in Fig. m. 

In addition to V2 — (cos(2(^)) (where the average is 
taken with the particle momentum distribution) we also 
show the p|.-weighted elliptic flow (denoted as a in S) 
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Fig. m shows that for pions V2 and W2^p2 differ by an over- 



all factor of about 2, but otherwise have the same impact 
parameter dependence. 
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FIG. 1. The 2"^* and 4"^ harmonic flow coeflicients for pions 
as functions of the number of participating nucleons. Also 
shown is the momentum-space anisotropy e^ (see text). 



This factor 2 is important: previously V2 and W2,p2 
have often been used synonymously and, based on 01- 
litrault's results [B|, one concluded that hydrodynamic 
calculations overpredict the elliptic flow at the SPS by 
about a factor of 2. Fig. shows that this is not the 
case: a correct comparison of the data with the calcu- 
lated V2 (not i'2,p2 !) shows good quantitative agreement. 
The data |1^ were obtained from Pb+Pb collisions at 
the SPS, with a cut on the collision centrality and on 
the particle px as given in the figure. Our calculation 
was done for 6=7 fm and, using Eq. (ra), the same px-cut 
as in the data was applied. The resulting values for V2 
at midrapidity are 2.9% for pions and 11.7% for pro- 
tons. The calculated pT-dependencc of V2 (not shown) 
also agrees with the data |£2[ , up to a normalization fac- 
tor which takes into account that we compute V2{pt) at 
midrapidity while the data [Q were obtained in 4 < y < 5 
where V2 is about a factor 3 smaller (see Fig. ||). 

The good agreement of the data (the shape of the pt- 
spectra as a function of b and the absolute values and 
PT-dependences of V2) with hydrodynamic calculations 
strongly suggests very early thermalization and pressure 
buildup in these collisions. In the calculation we can fol- 
low the time history of the elliptic flow: we found that 
W2.p2 for pions is nearly identical to 



{{Tx 



T, 



(5) 



yy/ 



if the spatial average 1 



is performed at the time when 



the fireball center freezes out. €p is the momentum-space 
analog of the spatial eccentricity ex defined above; it does 
not require knowledge of the particle spectra and can be 
evaluated also at other times from the solution of the 



hydrodynamic equations. As expected we find that e^ 
saturates as soon as the spatial anisotropy ex goes to zero. 
For Pb-|-Pb collisions at the SPS, ^ of the final elliptic 
flow is created while the fireball center is in a pure QGP 
phase, ^ of it is created in the mixed phase, and the final 
I is generated during the hadronic stage. This agrees 
with Sorge's conclusion WM that the elliptic flow at the 
SPS indeed probes the deconfining phase transition and 
the existence of a QGP phase. 
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FIG. 2. Elliptic flow V2 for pions and protons, as a function 
of rapidity, as measured by NA49 in 158 A GeV Pb-|-Pb colli- 
sions U^- The dashed lines are to guide the eye. The circled 
crosses at midrapidity show our hydrodynamic results, with 
the same cuts in b and pt as the data. 



However, is it also sensitive to the existence of a phase 
transition? To answer this question we recalculated V2 
and V4 with EOS I and EOS H, readjusting the initial 
conditions to the measured h~ and p ~ p spectra from 
central Pb+Pb collisions H. (While for EOS H an ac- 
ceptable fit is possible, the fit for EOS I is quite bad, as 
found before by several other authors.) Whereas EOS I 
(which can already be excluded from the 5 = spectra) 
gives about 30-40 % larger values for V2 , the elliptic flow 
developed by EOS H is quite similar to that of EOS Q. 
V4 is about 60% larger with EOS H than with EOS Q. 
The time history of e^ reveals that the softening of the 
EOS near the phase transition delays the buildup of ellip- 
tic flow by about 1.5-2 fm/c but that, at this particular 
beam energy, in the end it reaches the same value. The 
mechanism is the same as discussed in the context of van 
Hove's plateau: the phase transition weakens the elliptic 
flow, but since the system also spends more time in the 
transition region, its net effect on V2 is much less than 
naively expected. 

Given the apparent insensitivity of elliptic flow to the 
phase transition at a fixed beam energy, one may still 
hope for distinctive features in the excitation function 
of anisotropic flow [pO[. In Refs. |13^,H it was suggested 
that, for initial conditions around the "softest point", the 



response ^ r^ 2 ^ should develop a plateau-like struc- 
ture: as the beam energy or the collision centrality is 
increased and the initial energy density in the fireball 
center rises from subcritical to supercritical values, one 
should see a weaker elliptic flow response in the tran- 
sition region. To check this expectation we calculated 
excitation functions for U2.4 for Pb+Pb collisions at fixed 
impact parameter 6 = 7 f m (Fig. 0) . Since we only varied 
Co (the initial central energy density in 6==0 collisions), 
leaving the other model parameters (tq, eo/no, and Tdoc) 
unchanged, the curves in Fig. ^ correspond to constant 
initial eccentricity ex — 0.25. Hence they give directly 
the -^/s-dependence (parametrized via dN/dy\y^^) of the 
response of the anisotropic flow to ex ■ 
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FIG. 3. Hydrodynamic excitation functions of tlie 2" and 
4"^ harmonic flow coefficients, 112 and U4, for pions from A + A 
collisions {A ~ 200) at impact parameter b—7 fm. The verti- 
cal daslied lines indicate the produced total pion multiplicity 
densities at midrapidity for 11.6 and 158 A GeV beam energy 
(upper ends of the AGS and SPS ranges). The dash-dotted 
vertical line indicates the threshold above which, at b=7 fm, 
the fireball center is initially in a pure QGP phase. 



Fig. |3| covers the range 1 < Cq < 25 GeV/im^. Above 
SPS energies, V2 and v^ are seen to approach constant 
asymptotic values. This reflects the saturation of the 
anisotropic flow before freeze-out and its character as 
an early stage signature. V2 saturates before W4: first 
the dominant elliptic spatial deformation disappears, the 
smaller higher-order deformations are washed out later. 
At lower energies, V2 and W4 drop due to the decreas- 
ing fireball lifetime: the initial central temperature ap- 
proaches Tdcc7 causing decoupling before the anisotropic 
flow could fully develop. While this is a generic feature 
we must caution that quantitatively our results become 
unreliable in this region: We kept T^cc fixed although 
at lower beam energies freeze-out is known to occur at 
lower temperatures [ pl| ; this gives more time for flow 
buildup, leading to larger values V2 > 0. On the other 
hand, below 1-2 A GeV/c the elliptic flow at ycm builds 



up before the spectator nucleons have moved out of the 
way (as we assume); this "inertial confinement" causes 
the elliptic flow to develop perpendicular to the reaction 
plane ( "squeeze-out" 0] , f 2 < 0) instead of in-plane as in 
our calculations. Experimentally this sign change of V2 
occurs near i^bcam — ^A GeV [p2| . Neither of these two 
phenomena is, however, directly related to the existence 
of a phase transition, and above AGS energies our results 
are not affected by them. 

Fig. shows that the expected weakening of the 
anisotropic flow due to the phase transition sets in be- 
tween AGS and SPS energies: comparing the curves for 
the hadron resonance gas equation of state (EOS H, 
dashed) with those for EOS Q (solid), one sees that above 
SPS energies both V2 and v^ are reduced for EOS Q. At 
high energies the anisotropic flow cocfhcicnts thus show 
a qualitatively similar dependence on the EOS as the ra- 
dial flow [|2| . The reduction for V2 is relatively small (for 
pions V2 drops from 6 % to 5 %) ; W4 is reduced by about 
40% but, since Vi is so small, this is harder to measure. 
The effect on V2 is also more characteristic, due to the 
"bump" in V2 between AGS and SPS energies resulting 
from the interplay of the various effects discussed above; 
for U4 the softening of the EOS cannot break the basic 
monotony of the excitation function. 

Within hydrodynamics, the basic phase transition sig- 
nature in the anisotropic flow is thus a hump (and not the 
conjectured |13 33] plateau and second rise) in the exci- 
tation function of V2- Instead of changing -y/s, one can 
also vary the impact parameter at fixed ^/s to change the 
produced multiplicity density <iA^7r/dy|jj=ycm- This pro- 
vides for an alternate way to study the curves shown in 
Fig. 0, albeit not at fixed spatial ellipticity ex- We found 
p2[ that in this case the elliptic flow response ^ « 2^ 
develops a similar bump p4| . Unfortunately, for 158 A 
GeV Pb-fPb collisions it lies at the limit of the acces- 
sible impact parameter range (near 6=11 fm) where the 
hydrodynamic approach is expected to break down [B5|; 
for the planned SPS run at 40 A GeV, however, the bump 
should be clearly visible in semicentral Pb-t-Pb collisions. 

As a last point we discuss why the seemingly so dra- 
matic phenomenon of the "cracked nut" , recently advo- 
cated by Teaney and Shuryak |15| ] as a hydrodynamic 
signature for the existence of a QGP phase transition, 
doesn't leave stronger traces in the anisotropic flow pat- 
tern. These authors argued that at high energies (RHIC 
or LHC) a soft region in the EOS leads to the develop- 
ment of a "shell" at the edge of the almond-like initial 
fireball which is then cracked by the high pressure inside, 
with two separating half-shells expanding into the reac- 
tion plane. While we confirm their numerical results p5|, 
we tend to interpret them more cautiously. To illustrate 
our point of view we show in Fig. the freeze-out surface 
t{x, y) for a Pb+Pb collision at 6=8 fm, initiated with 
a central temperature To=870 MeV at to=0.2 fm/c and 



yielding a pion midrapidity density dN^r / dy\y=y^^f^5S0 
(corresponding to dN^^ / dy\y=y^^f»1600 for central colli- 
sions). Note its mushroom-like structure: While at the 
SPS hydrodynamics predicts freeze-out surfaces which 
shrink with time, the present surface features dramatic 
transverse growth pq | before freezing out nearly in- 
stantaneously after 13 fm/c. One thus expects a very 
small emission time duration signal in two-particle Bose- 
Einstein correlations ^^, in spite of the phase transition. 
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FIG. 4. Freeze-out hypersurface T{x,y,z—0) for b=8 fm 
Pb+Pb collisions, for an initial temperature ro=1.14 GeV 
{dNn / dy\y=y^^ ~ 1600) in central collisions. Note the dra- 
matic transverse growths, followed by sudden freeze-out. 



Already before freeze-out the initial elliptic spatial de- 
formation has vanished. The ripple on the top of the 
mushroom near its outer edge in x-direction is the "nut 
shell" [|l5|: in a cut through the surface at r « 13 fm/c 
it shows up as a crescent-shaped half shell at a; ~ 7 fm. 
However, a mere 0.5 fm/c later, the matter in this shell 
has frozen out, too. For EOS H one finds a very simi- 
lar mushroom, but without the ripple at the edge. Since 
there is no qualitative difference in the momentum-space 
structure of the "shell" compared to the rest of the mat- 
ter, this explains why it is impossible to uniquely identify 
this particular structure by an anisotropic flow analysis. 
As suggested in ||l5[, two-particle correlations may be 
more promising, but require extensive studies. 

We thank P.V. Ruuskanen, M. Kataja, R. Venugopalan 
and P. Huovinen for allowing us to use their hydrody- 
namics code for central collisions and to modify it for 
the study of non-central collisions. We are indebted to 
H. Dobler, H. Heiselberg, A. Poskanzer and S. Voloshin 
for constructive remarks and acknowledge fruitful dis- 
cussions with E. Shuryak and D. Teaney. This work was 
supported by BMBF, DFG and GSI. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 

[10] 

[11] 
[12] 

[13] 
[14] 
[15] 
[16] 
[17] 



[18] 



[20] 
[21] 
[22] 
[23] 
[24] 
[25] 
[26] 
[27] 
[28] 
[29] 

[30] 

[31] 
[32] 

[33] 
[34] 



[35] 



[36] 

[37] 



J. Barrette et al., Phys. Rev. Lett. 73 (1994) 2532. 

H. Appelshauser et al., Phys. Rev. Lett. 80 (1998) 4136. 

M.M. Aggarwal et al., Phys. Rev. Lett. 81 (1998) 4087. 

F. Ceretto et al., Nucl. Phys. A 638 (1998) 467c. 

J.Y. OUitrault, Nucl. Phys. A 638 (1998) 195c. 

H. Sorge, Phys. Lett. B 402 (1997) 251. 

H. Stocker et al., Phys. Rev. C 25 (1982) 1873. 

J.Y. OUitrault, Phys. Rev. D 46 (1992) 229. 

CM. Hung, E.V. Shuryak, Phys. Rev. Lett. 75 (1995) 

4003. 

E.V. Shuryak, O.V. Zhirov, Sov. J. Nucl. Phys. 28 (1978) 

247; Phys. Lett. B 89 (1980) 253. 

L. van Hove, Phys. Lett. B 118 (1982) 138. 

M. Kataja et al., Phys. Rev. D 34 (1986) 2755. 

H. Heiselberg, A.-M. Levy, Phys. Rev. C 59 (1999) 2716. 

H. Sorge, Phys. Rev. Lett 82 (1999) 2048. 

D. Teaney, E.V. Shuryak, ^cl-th/9904o"oe . 

S.A. Voloshin, Y. Zhang, Z. Phys. C 70 (1996) 665. 

The data were taken from an update of those 

published in [pfj and can be found at the URL: 



_e^ 



th/9904018 



http://na49info.cern.ch/na49/Archives/Images/Publica- 

tions/Phys.Rev.Lett.80:4136-4140t l998 

J.D. Bjorken, Phys. Rev. D 27 (1983) 140. 

H. Dobler, J. SoUfrank, U. Heinz, nucl- 
Phys. Lett. B, in press; E. Schnedermann, 
J. SoUfrank, U. Heinz, Phys. Rev. C 48 (1993) 2462. 
Particle Data Group, Eur. Phys. J. C 3 (1998) 1. 
J. SoUfrank et al., Phys. Rev. C 55 (1997) 392. 
P.F. Kolb, J. SoUfrank, U. Heinz, in preparation. 
F. Cooper, G. Frye, Phys. Rev. D 10 (1974) 186. 

Appelshauser et al., Phys. Rev. Lett. 82 (1999) 2471. 
Kampfer et al., J. Phys. G 23 (1997) 2001. 
Appelshauser et al., Eur. Phys. J. C 2 (1998) 661. 
Tomasik, U.A. Wiedemann, U. Heinz, in preparation. 
SoUfrank, P. Koch, U. Heinz, Z. Phys. C 52 (1991) 593. 
Braun-Munzinger, J. Stachel, Nucl. Phys. A 638 (1998) 



3c. 

P. Danielewicz et al., Phys. Rev. Lett. 81 (1998) 2438. 
N. Herrmann, Nucl. Phys. A 610 (199 6) 49c. 
C. Pinkenburg et al. |nucl-ex/9903010| , Phys. Rev. Lett., 
in press. 

E.V. Shuryak, Quark Matter '99, Nucl. Phys. A, in press. 
Our results suggest that the plateau in ^ found in jlj] 
must be a non-equtlibrium effect. Sorge himself stressed 
|l4{ that within his approach non-equilibrium mecha- 
nisms are responsible for the saturation of V2 as a function 
of time. An experimental study of the excitation func- 
tions (Fig. m) will thus provide valuable evidence for the 
degree of thermalization reached in the collision, through 
both their form and magnitude. 

Preliminary NA49 data on the 6-dependence of V2 do not 
show a bump (A. Poskanzer, Quark Matter '99, Nucl. 
Phys. A, in press). 

M. Kataja et al., Z. P hys. C 55 (1992) 153. 
U. Heinz, B.V. Jacak, |nucl-th/9902020| , Ann. Rev. Nucl. 
Part. Sci. 49 (1999), in press. 



